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Abstract 

We present a scaling investigation of some correlation functions in 0(a) improved 
quenched lattice QCD. In particular, as one observable the renormalized PCAC quark 
mass is considered. Others are constructed such that they become the vector meson 
mass and the pseudoscalar meson decay constant when the volume is large. For the 
present discussion we remain in intermediate volume, (0.75 3 x 1.5) fm 4 with Schrodinger 
functional boundary conditions. By fixing the 'pion mass' and the spatial lattice size 
in units of the hadronic scale ro, we simulated four lattices with resolutions ranging 
from 0.1 fm to 0.05 fm and performed the extrapolation to the continuum limit. The 
maximal scaling violation found in the improved theory is a ~ 6 % effect at a ~ 0.1 fm. 
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1 Introduction 



One of the major drawbacks in the standard formulation of lattice QCD, induced by 
the Wilson-Dirac operator violating chiral symmetry at the scale of the cutoff 
reflects among others in the fact that the quark masses are not protected from additive 
renormalizations and that the leading lattice effects in physical matrix elements and 
amplitudes are proportional to the lattice spacing a rather than being of 0(a 2 ). 

Nowadays, a systematic approach based on the Symanzik improvement programme 
H,[|] is well established to permit a removal of these discretization errors of 0(a) in 
lattice QCD with confidence. It has been elaborated for on-shell quantities in refs. (5]|| 
and adds appropriate higher-dimensional operators to the action and fields in order to 
compensate for any correction terms at 0(a). 

Within this framework a mostly non-perturbative 0(a) improvement of the ac- 
tion and the quark currents as well as their renormalization has been achieved in the 
quenched case, where quark loops are ignored. The basic idea for the practical im- 
plementation of on-shell improvement to enable a numerical computation of these im- 
provement coefficients — pioneered by the ALPHA collaboration — is to exploit chiral 
symmetry restoration and certain chiral Ward identities from Euclidean current algebra 
relations on the lattice at 0(a) [0j8|. Most prominently, demanding the PC AC relation 
d^A^ = 2mP a between the isovector axial current and pseudoscalar density P a to 
hold as a renormalized operator identity, imposes a sensible improvement condition in 
this respect. As a result, the values of the improvement coefficients c sw , ca and cy, 
which are required to completely eliminate the corresponding 0(a) effects, are known 
for P = 6/<? 2 > 6.0, go being the bare gauge coupling p|,10|,ll|. 

Therefore, one is now interested in the quality of scaling behaviour and the size of 
its possible violation. During the last two years it was reported in this context 12 . 13 , 14] 
that at a ~ 0.1 fm the residual 0(a 2 ) lattice artifacts may be still fairly large e.g. for 
the kaon decay constant fK r o (~ 10 %), while they are very small already for other 



quantities like the rho meson mass m p /^/a (~ 2 %) [14,15]. Thus, restricting to an inter- 
mediate physical volume, the present scaling tests are intended to examine the impact of 
full 0(a) improvement in quenched lattice QCD thoroughly and with high accuracy for 
different observables. Some of these are designed to coincide with phenomenologically 
relevant observables in the limit of large physical volumes. A preliminary status of this 
work has been briefly surveyed recently in ref. [16]. 

We wish to point out that the present investigation must be distinguished from sim- 
ilar studies like [13,17,18] for an important reason. Namely, since we address directly 
the scaling behaviour of correlation functions in finite volume, we do not have to rely on 
the asymptotic behaviour of (ratios of) usual timeslice correlation functions to extract, 
for instance, hadron masses or decay constants; these more conventional techniques of- 
ten reveal to be genuinely affected by systematic errors difficult to control. Actually, 
an avoidance of such intrinsic uncertainties is supplied in part by the Schrodinger func- 
tional: its finite- volume fermionic correlation functions are scaling quantities at any fixed 
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distance in time, if the renormalization factors of the quark fields at the boundaries are 
properly divided out. Beyond that, they decay very slowly for small time separations, 
allow to gain a good numerical precision and hence offer the appealing possibility to 
probe the theory for 0(a) improvement without unwanted additional sources of errors. 
Finally for sufficiently large volumes the correlators can be shown to embody standard 



hadronic masses and matrix elements [19|. These aspects provide the real advantages of 
our method. The prize to pay for it is, however, that within the Schrodinger functional 
formulation there exist two further improvement coefficients multiplying the boundary 
counterterms. These are perturbatively known only and as a consequence, one can in 
principle not exclude that the theory is to some extent still contaminated with uncan- 
celled 0(a) contributions. That this fear can be dismissed in practice has recently been 
demonstrated for the renormalization group invariant running quark mass pQfl , but we 
will make sure of it in the present context too. 

The paper is organized as follows. In section |2| we introduce the fermionic corre- 
lation functions under study and sketch how spectral observables can be constructed 
from them in the Schrodinger functional scheme. After a short account on the numer- 
ical simulations, section || contains a detailed description of the scaling test and the 
careful evaluation of the data. Section |] gives the results. Here, also the question will 
be answered whether an improvement condition, chosen to fix a certain improvement 
coefficient non-perturbatively, is unambiguous in the sense that it automatically im- 
plies appreciably small scaling violations of 0(a 2 ) in other quantities not related to this 
specific condition. We conclude with a discussion in section || 

2 Correlation functions and hadronic observables 

The basic framework for our lattice setup is the QCD Schrodinger functional (SF), whose 
concepts and characteristic features have been published in much detail in refs. |21,22, 



23]; consult also [^J,^] for comprehensive overviews on the subject. 

It is defined as the partition function of QCD in a cylindrically shaped space-time 
manifold of extension L 3 x T with periodic boundary conditions in the space directions 
and (in general inhomogeneous) Dirichlet boundary conditions at times xq = and xq = 
T. This means in the case of the gluons to require the spatial gauge field components 
at the boundaries to satisfy Ak(x)\ xo= o = C/^x) and Ak(x)\ X0= T = C' k (x), where Ck 
and C' k are some prescribed smooth classical (chromo-electric) gauge potentials, and a 
similar assignment is imposed on the quark fields as well. One of the bold advantages 
of such a choice is that it ensures frequency gaps for the gluon and quark fields, and 
thereby numerical simulations at vanishing quark masses become tractable. As in most 
of the other applications of the QCD SF, we assume the special choice of homogeneous 
boundary conditions from now on: = C' k = for the spatial components of the gauge 
potentials and vanishing fermion boundary fields. 



2 



2.1 Correlation functions in the Schrodinger functional 

Although the definitions of fermionic correlation functions within the SF already ap- 
peared in the literature 1^,11], let us recall them here and collect the essential 
properties and formulae in order to make the paper self-contained. 

If C and C denote 'boundary quark and antiquark fields' at Euclidean time xo = 
and primed symbols the corresponding objects at Xq = T [j?fl, one builds up the boundary 
field products 

O a = « 6 E C(y)75 \ r a CM , O fa = a 6 £ C'(u)75 \ r a C(v) (2.1) 
y.z u>v 

and analogously, 

Ql = « 6 E C(y)Tfc 2 r a CM , <# = ^ £ CVh* | r a C'(v) , (2.2) 

y,z u,v 

where T a , a = 1,2,3, are the Pauli-matrices acting on the first two flavour components 
of the quark fields if). In the operator language of quantum field theory they create 
initial and final quark-antiquark states with zero momenta, respectively, and transform 
according to the vector representation of the exact isospin symmetry. For the axial- 
vector current and the pseudoscalar density P a we use the local expressions 

A%{x) = ?(*) 7m 75 § r a i>{x) , P a (x) = V^) 75 | r a i>{x) , (2.3) 

while the vector current and the anti-symmetric tensor field read 

V£(x) = ?(x) 7m \ r a i>(x) , T^(x) = i^(x)o^ \ r a ^(x) . (2.4) 

Now we consider correlation functions on the lattice in the SF. By inserting the preceding 
densities at some inner point x of the SF cylinder (with support on the hypersurface 
at xq) between the appropriate external quark-antiquark states, one introduces the 
expectation values 



/ A (xd) = "« 6 E^ (^WC(y)75h a (W) = -\{A a {x)O a ) (2.5) 
/p(x ) = -a"Y,\ <^ a (^)C(y)75|r a CM) = -\ (P a "(x) O a ) , (2.6) 
analogously, 



3 x v-/^w//-2 -v"// 3 
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and the boundary-boundary correlation function 



ky(x ) = -I {V k a (x)Qi) (2.7) 
^t(xo) = -I (TZ (x)Q a k ) , (2.8) 
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Gauge invariant correlators of this type have already been used to study the conservation 
of currents on the lattice and to deduce suitable improvement and normalization condi- 
tions in lattice QCD in order to calculate the corresponding coefficients and constants 
non-perturbatively by numerical simulations ||.1C.11]. One of them, /i, will be utilized 
later to cancel the multiplicative renormalization of the boundary quark fields 
Besides on the SF characteristic kinematical variables 0, the correlation functions de- 
pend on the bare parameters go, mo and the improvement coefficient c sw = c sw (go) in 
the fermionic part of the lattice action, but not on the spatial coordinates of x owing 
to translation invariance. There is also a dependence on the improvement coefficients 
ct and c"t, which account for specific boundary 0(a) counterterms arising in the SF 
approach Q. 

After contracting the quark fields, all the bare and unimproved correlation func- 
tions have the general structure hx(xo) °c (Tr {H + (x)TxH (x)}), h = f,k, where the 
respective insertions are Tx £ {— 70, 1, 7fc, — o ko}-, X=A,P,V,T; the trace extends over 
colour, Dirac and (in principle as well) flavour indices, and the matrix H is the quark 
propagator from the boundary at xq = to the point x in the interior of the space-time 
volume ||. In the quenched approximation the expectation values are understood to 
be taken as path integral averages in the pure gauge theory. 

We should mention that also the primed correlation functions h' x , which are con- 
nected to hx through a time reflection and vice versa, become relevant. Those are in 
the case of eqs. J2.5|) and ( [2.6D 



f A (T - x ) = +1 (A a Q (x) 0' a ) , /p(T - so) = -~ (P a (x) O la ) , (2.10) 

and similar relations apply to ky,ky and kT,k' T - Obviously, in h' x the currents and 
densities are probed by the boundary quark fields at xo = T instead, and the argument 
T — xo indicates that they fall off with this distance. For vanishing gauge fields C and C 
at the boundaries, our correlation functions possess the useful time reflection invariance 
hx(xo) = h' x (xo)- This allows to sum them up accordingly, and averaging over the 
spatial components helps to reduce the statistical noise in the Monte Carlo simulation 
further. 

On-shell improvement at 0(a) for the axial and vector currents is achieved by 
adding the derivatives of the pseudoscalar density and the tensor current as the suitable 
0(a) counterterms, 

(Aj^(x) = Al{x) + ac A d^P a (x) (2.11) 

(Vi) a ^x) = y;(x) + ac v a^(x), (2.12) 

where the improvement coefficients c\ and cy are determined by the demand to can- 
cel the 0(a) errors in lattice Ward identities, emerging from a mixing with higher- 
dimensional operators with the same quantum numbers [10,11 1. Then the corresponding 
improved fermionic correlation functions are given by: 

A( x o) = fA(xo)+ac A d fp(xo) (2.13) 
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ky(x ) = ky(x ) + acvdohrixo) ■ (2.14) 

The lattice derivative = ^ (d^ + <9*) is the symmetrized combination of the usual 
forward and backward difference operators and <9*, acting as 

dj{x) = f{ - x+a ®- f{x \ dim = f{x) ~ /(x = qA) . 

a M a 

Herewith we are already in the position to write down the unrenormalized PCAC quark 
mass as a function of the timeslice location a < xq <T — a: 

= d f A (x ) + ac A d*d f P (x ) 
2/p(x ) 

For the properly chosen value of c A = ca(<?o) at given gauge coupling g and any 
hopping parameter k it is defined by obeying the PCAC relation (for two degenerate 
quark flavours) up to cutoff effects of 0(a 2 ), 

d II [(A l );(x)]=2mP a (x) + 0(a 2 ), (2.16) 

which more rigorously must be looked at as a renormalized operator identity 

d, ((A R )«(x)0) = 2m R (P^(x)O) + 0(a 2 ) 

in terms of some arbitrary renormalized on-shell 0(a) improved field O localized in a 
region not containing x. 

2.2 Construction of renormalized observables 

Now we want to introduce the renormalized correlation functions and design the scaling 
combinations of them, which will be studied numerically in the next section. 

The QCD SF serves as a particular intermediate finite-volume renormalization 
scheme which is, however, not necessarily related to a special regularization j7,24,25]. 



Here, the SF is employed as a mass-independent renormalization scheme, while the ra- 
tio T/L is assumed to be kept fixed to a certain value. The freedom in choosing the 
boundary fields C and C (as well as the boundary conditions on the quark fields spec- 
ified by angles 9^, cf. subsection [O]) different from zero are left for other applications, 
see e.g. ||,2C,26]. Moreover, the SF respects 0(a) improvement after adding the 0(a) 



counterterms oc c t , ct so that by attaching additive and multiplicative renormalization 
constants, the quantities 

(A R )-(x) = ZAil + bAam^A^x) (2.17) 
(Vr)«(x) = Z v (l + 6 v am q )(y I )^(x) (2.18) 
P&(x) = Zp(l + b P am q )P a (x) (2.19) 
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induce the renormalized and improved correlation functions 



/£(x ) = Z A (l + & A am q )Z 2 (l + 6 c am q ) 2 / A (:ro) (2.20) 

k§(x ) = Z w (l + byam q )Z^(l + b ( am q ) 2 k l y(xo) (2.21) 

f$(xo) = Zp(l + 6 P am q )Z 2 (l + 6 c am q ) 2 /p(2;o) (2.22) 
-r v^urrtq^ 



/ X R = Z f 4 (l + Um q ) 4 /i • (2.23) 



The constants and 6^ (the former being scale dependent J7j) have to be attributed to 
the boundary values of the quark and antiquark fields appearing in the products (2.1) 



and (^^) in the renormalized theory. In a mass-independent renormalization scheme 
the underlying 0(a) counterterm enters as 

Cr(x) = Z C {1 + fr ( am q )C(x) 

and similarly for the antiquark field £, giving 0^ = Z 2 (l + b^am q ) 2 O a for instance. 

Because the renormalization does not distinguish between different flavours in a 
mass-independent scheme, the knowledge of the Z-factors suffices to link the lattice 
theory at finite cutoff to the renormalized continuum theory. Here, all normalization 
conditions to fix and determine the renormalization constants Zx and 6x, X=A,V,P, 
were imposed on appropriate matrix elements at zero quark mass, which is safe within 
the SF scheme as the finite lattice extent L provides the before-mentioned natural 
infrared cutoff for the theory |}],|lO| . Because of the zero quark mass condition they are 
functions of go only and not of the subtracted bare quark mass, which equals 

am q = am - am c (g ) = ~ ( — ) (2.24) 

and vanishes along a critical line thq = m c (go) in the plane of bare parameters, implicitly 
defining the critical hopping parameter k c . Any remaining corrections of 0(am q ) are 
supposed to be cancelled by adjusting the bx alone. 

Now we can pass to the set of observables we have constructed for the present 
study. We start with the renormalized PCAC (current) quark mass in the SF scheme, 
which in view of eqs. ( |2.15| ), ( |2.17| ) and ( |2.19| ) may be defined as 

m = ^™(f) (2.25) 

by multiplying the proper renormalization constants Z& and Zp(L), the latter assumed 
to be taken at some renormalization scale fi = 1/L [|20| , |27| . Strictly speaking, the ratio 
of the additive renormalization factors 1 + b^am^ and 1 + bpam q would have been to be 



accounted for as well, but it turns out perturbatively J28|] and non-perturbatively |29[] 
that (6a — bp)am q is numerically quite small at the interesting values of the bare gauge 
coupling and the hopping parameter. Hence we neglect it, here. 
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Beyond that, we compose the following (time dependent) combinations of renor- 
malized and improved fermionic correlation functions. Firstly, the logarithmic time 
derivatives 

doff(x ) d f P (x ) 
m(X0) = If&oT = ~MxT (2 ' 26) 

Mv ( Xo ) = g o fc v( x o) = goM^o) (2.27) 
ky{xo) k v (x ) 

of the respective SF correlation functions in the pseudoscalar and vector meson channel; 
they deviate from the ordinary definition of effective masses (~ do In hx) by terms of 
0(a 2 ). Secondly, we will consider the ratios 

/ ^ r ^a(^o) 
^ps(^o) = Cp S 



'/f 

C PS ZA(1 + h 2^M^l , Cpg= 2 (2.28) 



w(^o) = cv- 



7f 

Since through the division by ^/ f^ - the renormalization factors of the boundary quark 
fields, Zq(1 + b^arriq), drop out, it is ensured that tjps an d f]y ( & t fixed argument 
xo) exhibit scaling and have a well-defined continuum limit. We note in passing that 
alternatively the correlation function f\ also could have been used instead of fp to 
define a 'local mass' in the pseudoscalar channel. But as f\ amounts to somewhat 
larger statistical errors in the time derivatives, we here preferred fp for the purpose of 
Vps(x ). 

In the present scaling test we fix a definite temporal separation from the boundaries, 
xq = T/2. Thus one arrives at the objects 

pseudoscalar channel : [ips (^) rjps (^) (2.30) 

vector channel : fly (-0 rjy ( -^) ■ (2-31) 

Their continuum limits are expected to be approached like /ix(^-) + 0(a 2 ) and ^x(f-) + 
0(a 2 ), X=PS,V, in the 0(a) improved theory. The choice xq = T/2 is motivated by 
the fact that in the SF scheme cutoff effects are generically larger when sitting closer 
to the boundaries. It is possible to construct other quantities in a similar way, but we 
consider the foregoing ones as reasonably representative. 
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Let us emphasize, however, a final point. Adopting the quantum mechanical rep- 



resentation of the field operators associated to (2.3) and ( |2.4| ), it can be shown in the 
transfer matrix formalism that asymptotically for large Euclidean times the quantities 
in eqs. (|2.30|) and ( |2.31 ) become the pseudoscalar and vector meson masses (^tps,Mv) 



as well as the pseudoscalar decay constant (??ps)- For instance, the proportionality con- 
stant Cps in ( [2,28 ) is such that this ratio turns, as the temporal lattice extent goes 



to infinity (xo,T — xq —* oo), into a familiar matrix element, which complies with the 
standard definition of the pion decay constants in continuum QCD: 

Z A ((#(x) 7 o75 5 r a V(x)|vr a (0)) = m*/* ■ 

More formally, the observables O just introduced should be regarded as functions 
0(T/L,xo/L, L/ro,a/ro), and in the spirit of the above a physically meaningful situa- 
tion is realized if, as the spatial volume (L/a) 3 tends to infinity, xq 3> ro and T—xq r$ 
are valid for some typical hadronic radius of tq ~ 0.5 fm. Masses and matrix elements 



of interest in hadron phenomenology can then be extracted [19]. 
3 Scaling tests 

Before going into the details of the investigation, one has to keep in mind that here 
the actually taken lattice volumes in physical units are only of intermediate magnitude. 
Therefore, any of the following results are prevented from resembling the large volume 
limit, where our observables were argued to receive a really physical meaning. Instead of 
this, most emphasis is on the scaling properties of the theory, which should not depend 
on the specific choice of the lattice size and the SF characteristic boundary conditions. 
Nevertheless, we refer from now on to the (first three of the) quantities in eqs. ( |2.30D 
and ( 2.31| ) as the 'pion mass', its decay constant and the 'rho meson mass' and assign 
the common symbols m n , and m p to them in the obvious manner. We also denote 
fjy = ??v(^") m the vector meson channel. 

The advantage of working at finite quark mass in a direct test of improvement 
should be stressed explicitly. Namely, this avoids any extrapolations and evades the 
potential problems with so-called exceptional configurations one runs into, when the 
parameter region of zero quark mass is attempted to be reached p|,^Q]. 

3.1 Monte Carlo simulation 

The cost of a quenched QCD simulation is always governed by the computation of 
fermion propagators required for the correlation functions to be covered. This involves 
the action of the Wilson-Dirac operator on quark fields tp(x), which for the 0(a) im- 
proved theory in the framework of the SF is conveniently decomposed as D + 5D + tuq 
with 

(D + m )^(x) = ^-M^(x), K = 1 (3.1) 
2k 8 + Zitiq 
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Mtp(x) = ij){x) - ^XfJJ(x,fj,)(l - Jfji/jjx + ay) 

+X* ti U(x - ay, y)+(l + 7m )^(x - a/i)} , (3.2) 



where, as usual, the fermionic degrees of freedom ip(x) live on the sites x of the lattice, 
and U(x,y) denotes the SU(3)-valued gauge links in lattice direction x + ay, y = 
0, . . . ,3. The factors = Q ta9 i j -/ L j Q = fj, give rise to a modified covariant derivative 
equivalent to demanding spatial periodicity of the quark field up to a phase e tdk ; for 
our purposes 8k, k = 1,2,3, was set to zero throughout.^ The local 0(a) counterterm 
5D = 5D V + 5Db now consists of two contributions, namely the Sheikholeslami-Wohlert 



clover term [31| 



i 

SD v ij}(x) = c sw - aa /JiV F IJiV (x)'il}(x) (3.3) 

and a term 

5D h ^{x) = (c t -l){5 Xo/aA liP(x)-U(x-a6,0) + P + iP(x-a6)\ 

+S xo/ajT ^(x)-U(x,0)P^(x + a6))} (3.4) 

with 

P + ^{x) =0, P^(x) m = 0, P± = i(l± 7o ), 



xo=0 



x =T 



which is specific for the SF type of boundary conditions in our setup. The improvement 
coefficient ct and a further one, ct which enters the calculation as well but is independent 
of the local composite operators containing quark fields and thus not written down here, 
are set to their one- loop perturbative values .0 

Since the technicalities of the Monte Carlo simulations are identical to those already 
detailed in ref. || , it is not necessary to repeat them here in full length. Our data were 
taken on the APE-100 massively parallel computers with 128 - 512 nodes at DESY 
Zeuthen, whose topology also allows to simulate independent replica of the system at 
the same time in the case of smaller lattice volumes (sets A - C below). The gauge 
field ensembles were generated by a standard hybrid overrelaxation algorithm, where 
each iteration consists of one heatbath step followed by several microcanonical reflection 
steps (typically A^or = L/2a + 1), and the correlation functions have been evaluated by 
averaging over sequential gauge field configurations separated by 50 iterations. To solve 
the system of linear equations belonging to the boundary value problem of the Dirac 
operator within the measurements of the correlators, the BiCGStab algorithm with 
even-odd preconditioning was used as inverter. Finally, a single-elimination jackknife 
procedure was applied to estimate the statistical errors of all the secondary quantities, 

1 In general, 9^ can serve as an additional kinematical variable to formulate proper improvement 
and normalization conditions, see for instance [p|,po|. 

2 Recently, the coefficient has also been computed up to two-loop order of perturbation theory in 
the quenched case |B2| . 
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because the data stemming from the same configurations must be considered as strongly 
correlated. By dividing the full ensemble of measurements into bins we also checked for 
the statistical independence of our data samples. 



3.2 Method and numerical analysis 

For the analysis we use c sw , ex, Zx, X=A,V,P, and by non-perturbatively determined 
in [p|JTojJTT , 20 ] for (3 > 6.0, while 6a and bp are taken from one-loop perturbation theory 



0, | ,25]. If available, we always adopted the rational formulas for the former with overall, 

1. e. statistical and systematic, uncertainties stated in the references. As opposed to Z\ 
and Zy (and as already anticipated below eq. ( |2.25p in subsection 2.2), the normalization 



constant of the pseudoscalar density, Zp = Zp(L), acquires a scale dependence through 
its renormalization. The scale evolution of Zp, which due to its definition at the point 
of vanishing PC AC quark mass (but in the absence of exact chiral symmetry at finite 
a) is unique only up to 0(a 2 ) errors, has been recently computed non-perturbatively 
in p0[ . In the SF the appropriate renormalization scale is /j, = 1/L = l/2L max , with 
L max /ro = 0.718 and we take over the needed numbers for Zp from the last but one 
reference. In the case of cy and the critical hopping parameter k c in eq. ( |2.24D , where 
no closed expressions are recommended yet, we adapted the numbers at discrete values 
of the gauge coupling from refs. HP] to our simulated /3-values by linear interpolation. 

The strategy was then to keep a finite physical volume and the quark mass fixed by 
prescribing the geometry T/L = 2 and two further renormalization conditions, which 
we decided to chose as 

m n L = 2.0 and — = 1.49 (3.5) 
ro 

for the 'pseudoscalar meson (pion) mass' and the spatial lattice size, respectively. The 
first condition on can be, at least approximately^], satisfied by a careful tuning 
of the hopping parameter (a 2 m 2 ~ k). Here the reference scale is expressed by the 
hadronic radius ro defined in []34j through the force between static quarks to yield the 
phenomenologically motivated value of ro — 0.5 fm. Using the latest results on the 
hadronic scale ro/a in ref. [33|, quoted there as 



In 1-1 = -1.6805 - 1.7139(/3 - 6) + 0.8155(/3 - Q) z - 0.6667(/3 - 6) d , (3.6) 

one can solve numerically for (3 after inserting a/ro = 1.49a/L to find the desired pairs 
(L/a,{3) in order to fulfill the second condition in eq. (|3.5| ) within errors. In practice, 
the particular value L/tq = 1.49 was determined by the initial simulations at [3 = 6.0 
on lattices with spatial size L/a = 8, and the larger lattices were adjusted thereafter. 
The simulation parameters and some results are compiled in table |l[ These settings 
give an intermediate volume of (0.75 3 x 1.5) fm 4 , and one moves on a line of constant 

3 In our actual simulations we have some mismatch in m n L between the individual points in param- 
eter space, which will be discussed later. 
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set 


L/a 


(3 


K 


-^meas 


L/r 


am n 


m n L 


A 


8 


6.0 


0.13458 


12800 


1.490(6) 


0.2505(11) 


2.004(9) 


B 


10 


6.14 


0.13538 


3840 


1.486(7) 


0.1945(14) 


1.946(14) 


C 


12 


6.26 


0.13546 


2560 


1.495(7) 


0.1709(13) 


2.050(16) 


D 


16 


6.48 


0.13541 


3000 


1.468(8) 


0.1244(9) 


1.991(15) 



Table 1: Simulation points and its statistics N meas denoting the number of gauge Held 
configurations, on which the fermionic correlation functions were computed. L/r$ and 
m v L are the quantities chosen to fix renormalization conditions for the LCP studied. 



physics (LCP) in bare lattice parameter space with lattice resolutions ranging from 0.1 
fm to 0.05 fm. 

As an important prerequisite for the reliability of the scaling test we had, of course, 
to estimate the dependence of our results on the SF specific (and solely perturbatively 
known) improvement coefficients of the boundary counterterms Ct and Ct, which with 
respect to full non-perturbative O(a) improvement represent the only imperfectly known 
input parameters for the simulation. Otherwise a complete suppression of errors linear 
in a would not be guaranteed, and continuum limit extrapolations with an 0(a 2 ) term 
as the dominant scaling violation is a priori not justified. To this end we verified by an 
artificial variation of the one-loop coefficients in the expansions 

mm 

c i-ioo P = !_ cosggg + OGri) (3.7) 
-i-ioop = 1 _ QM8g 2 + 0(fl 4 } (3 8) 

by a factor 2 for c\ loop and by a factor 10 for loop that at unchanged renormalization 
conditions (|3.5| ) their influence on the level of numerical precision in our data is small 
enough to be neglected: it typically came out to be below 1 % for af n , below 2 % for 
r/y, and nearly not visible for am and am p , at parameters corresponding to simulation 
point A (T/a = 16). It was sufficient to do these replacements for the lattice with the 
largest a, since the relative contribution of the boundaries to the field variables residing 
on the bulk of lattice points among a given gauge field configuration decreases with 
decreasing lattice spacing. Therefore, the leading scaling violations can be regarded as 
being purely 0(a 2 ) within our precision. 

In figure |l| we first illustrate for simulation points A and D the correlation functions 
/A; /p and some of the quantities deduced from them in the previous section in depen- 
dence of the time coordinate. The correlators reflect that a good signal remains also at 
larger distances in time. The PC AC quark mass ( |2.15 ) already exposes a plateau for the 
rather moderate temporal extensions of the lattice, whereas the 'local pion mass' ( p. 26 ). 
and even more the logarithmic derivative of the improved correlation of the axial current 
( [2.13 ), show significant contributions from higher intermediate states in small volumes. 
On physically large volumes, however, this pattern disappears: plateaux develop around 
xq = T/2, i.e. the lightest excitations, which coincide with the pseudoscalar meson mass 
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Figure 1: Correlation functions /a and fp (upper parts) and 0(a) improved local 
masses extracted from them: PCAC current quark mass via eq. ( |2J5| ) (middle part) 
and 'pseudoscalar meson masses' via eq. ( 2.26J and the same with fp substituted by f A 
(lower parts). The left hgure depicts simulation point A (L/a = 8) and the right one 
point D (L/a = 16), with T/L = 2 in both cases. The solid lines are only meant to 
guide the eye, and the statistical errors are smaller than the symbols. 



when denned either by f\ or by fp, govern the exponential decay of these functions. 
This will be explicitly demonstrated elsewhere ]19[| . An analogous statement holds for 
the vector meson mass via the correlation functions k\j and fcp- 

The expectation values in the simulation points of table [l] for the observables 
eqs. (ETlfil) , (HH) and djfip are now summarized in table |2[ Within all potential 
sources of errors to be incorporated in the analysis, i.e. the statistical one and those 
coming from Z\/Zp, Za, Zy, L/tq and m^L, the contributions caused by the uncer- 
tainties of the renormalization factors Zx ever dominate the combined errors quoted in 
the second parentheses in the table. Furthermore, any inherent small mismatch of the 
central values with the renormalization conditions (3.5) on m n L and L/tq in sets B, C 
and D was corrected by a conservative estimation of the slopes dO/d(am 1T ) and dO/d/3 
with O = 0(m n L, L/rQ, L/a) £ {am,am p ,af n ,rjy}, which enter the identities for the 
required partial derivatives 

d ° ^ 1 (3.9) 

c^jTiTr-L) L d(am 7T ) 

dO ^a dp dO^rp dO 
d{L/r ) L d(a/r ) 5/3 a d(L/a) ' 1 ' ' 

They were numerically extracted in linear approximation from several simulations in set 
A at neighbouring values of n and (5, where in the latter case they had to be combined 
with the derivative d(3/d(a/ro) to be read off from the parametrization eq. (|3.6|). The 
slopes obtained in this way were carried over to the finer lattices as well, since for 
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increasing lattice resolution (/? > 6.0) their corrections are of 0(a) and, with respect to 
the other sources of errors, can safely be ignored. 



set 


HI) 


am p 




Vy 


A 


0.0132(2) 


0.3438(10) 


0.1248(5)(13) 


0.3283(19)(25) 


B 


0.0083(1) 


0.2688(15) 


0.1013(7)(12) 


0.3387(36) (40) 


C 


0.0096(1) 


0.2337(14) 


0.0851(7)(11) 


0.3289(40) (43) 


D 


0.00651(6) 


0.1689(10) 


0.0650(5)(8) 


0.3561 (41) (44) 



Table 2: Analysis results for the observables under consideration in lattice units. The 
Erst error is only the statistical one, the second one (where given) includes in addition 
the uncertainties from the renormalization constants. 



4 Results 



We pass to the final results. After the procedure of matching the conditions charac- 
terizing the LCP under study, one finds the numbers collected in table || note that 



according to ( 2.25 ) the bare PC AC quark mass Tn(-n) translates through multiplication 



with Z\/Zp into the renormalized quantity m. Now we are prepared to perform ex- 



set 


mro 


m p r 


fnT0 


Vv 


A 


0.1069(50) 


1.846(13) 


0.6701(82) 


0.3283(46) 


B 


0.1029(35) 


1.839(15) 


0.684(11) 


0.3296(72) 


C 


0.1045(36) 


1.848(18) 


0.681(13) 


0.3360(95) 


D 


0.1103(35) 


1.860(18) 


0.699(12) 


0.3457(75) 



Table 3: Dimensionless results for the quantities of table^ with total errors after they 
have been corrected to fulfill exactly the simultaneous renormalization conditions (SI 



as described in the text (subsection 3.2). The errors on r^/a quoted in [33] have been 
taken into account as well. 



trapolations of these data to the continuum limit, assuming convergence with a rate 
proportional to a 2 . 

The fits are displayed in figures ^ - || and exemplify the scaling behaviour on the 
course from simulation points A to D. One evidently observes the leading corrections 
to the continuum to be compatible with 0(a 2 ). Moreover, it can be inferred from table || 
that the differences of the continuum limits from the values at f3 = 6.0 (a — 0.1 fm) are 
around or even below 5 % in the improved theory. These appear to be partly smaller 
than it was to be expected on basis of the experiences reported previously in refs. [13. 14] . 

At this point we have to add the remark that the results in the vector channel had to 
be revised compared to those listed in |ffq]. Due to some incorrect normalization of the 
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Figure 2: Scaling behaviour of the renormalized PC AC current quark mass in units 
of tq. An intermediate volume with SF boundary conditions is considered. The renor- 
malization constant Zp refers to a scale of L = 1.436 tq. Non-perturbative 0(a) im- 
provement allows to extrapolate linearly with (a/ro) 2 — > to the continuum limit. 




correlation functions ky and &t during an earlier data analysis, we erroneously observed 
a quite steep slope when carrying out the (a/ro) 2 — ^ fit for the ratio rjy defined in 
eq. ( 2.31| ). Looking at the final numbers now, it does no longer stand in contradiction 
to the findings in the pseudoscalar channel. By contrast, since the scaling violations of 
rjy are only slightly larger than for f n ro, we interpret this as a further compelling and 
satisfactory indication for the effectiveness of non-perturbative O(a) improvement. 

In ref. [14] the suspicion was raised that scaling looks even slightly better, if the 
perturbative estimates for the improvement coefficients ca and cy are used. This issue 
deserves some comments here. In order to address the sensitivity of the analysis to the 
improvement terms in the pseudoscalar and vector channel proportional to c\ and cy, 
we just evaluated our data with setting these coefficients arbitrarily to zero. Then, of 
course, the theory is only partially O(a) improved, and a residual contamination with 
uncancelled O(a) contributions has still to be expected. The surprising outcome is that 
upon omitting those terms, two quantities (/•7rn)j'7v) h ave somewhat smaller a-effects 
in total magnitude. At the same time, however, the renormalized PCAC quark mass 
mro gets much larger ones. Nevertheless, as seen from table ^j, the continuum limits 
of both data sets agree within errors. Such a result might suggest that the qualitative 
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Figure 3: Continuum limit extrapolation as in figure^ but for the 'rho meson mass' 
in intermediate volume. 



scaling behaviour of our quantities (apart from mro) is only marginally influenced by the 
definite choice of cx, X=A,V, still meeting the condition of a dominant a 2 -behaviour. 
On the other hand, we observed a tendency in the ca = cy = data points to disperse 
around the straight line fits to the continuum, which hints at some remnant of admixture 
of O(a) discretization errors; hence a scaling violating term oc o 2 in leading order rather 
seems to be ruled out in that case. This is particularly pronounced for mro, where 
the correction to the continuum limit grows distinctly (from 2 % to 17 % in table |j) 
if the improvement of the axial quark current is switched off. Opposed to that, in the 
fully improved case including the O(a) correction terms, the required continuum limit 



case 


mr m p r Ur n y 


fully improved 


0.1092(45) 1.856(20) 0.704(13) 0.3472(84) 
2 % < 1 % 5 % 6 % 


ca = and cy = 


0.1069(46) 1.856(20) 0.698(13) 0.3439(85) 
17 % < 1 % 1 % 2 % 



Table 4: Continuum limits and their percentage deviations from (3 = 6.0 (a ~ 0.1 fm). 
The lower numbers belong to a data evaluation, where the improvement coefficients ca 
and cy have been artificially taken to vanish. 
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Figure 4: The same as in figure ^ but for the 'pion decay constant'. 



extrapolations are generically not critical. In conclusion — and as a further convincing 
argument for the use of its non-perturbative values — this finally supports the physical 
insight that ca and cy are essentially relevant for chiral symmetry restoration at finite 
cutoff. 

5 Discussion and outlook 

In confirmation of similar investigations [[1^, 13 , 14, lq], O(a) improvement implies a sub- 
stantial reduction of scaling violations. Our numerical simulations of renormalized cor- 
relation functions in intermediate physical volume within the Schrodinger functional 
give clear evidence for an overall behaviour completely consistent with being linear in 
o 2 at a < 0.1 fm, for all quantities under consideration. Changing a by a factor two 
yields very stable fits and honest continuum limit extrapolations. Actually, the residual 
0(a 2 ) cutoff effects at a ~ 0.1 fm stay around {f^r^ and TJy) or significantly below 5 % 
(mro and m p ro). 

To quantify directly the influence of the improvement coefficients ca and cy on the 
scaling behaviour, we examined also the partially improved case ca = cy = 0. In this 
case f n ro and fjy show an even weaker dependence on the lattice spacing, but now one 
finds ~ 17 % lattice spacing effects at a ~ 0.1 fm in the renormalized current quark mass 
fnrQ. Additionally, as outlined in section |||, the functional form of the leading a-effects 
then appears no longer compatible with a 2 alone; furthermore, chiral Ward identities 
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Figure 5: Again as in figure^ but now for the observable n v , which is composed of a 
(renormalized and improved) correlation function involving the vector current. 



are badly violated for ca = cy = at 0(a) level [^,|TT|. Thus there is in general no 
choice for the improvement coefficients ca and cy, which diminishes the size of 0(a 2 ) 
lattice artifacts simultaneously for all relations and observables below a level of 5 %. 
However, one should not feel tempted to judge this fact as a kind of principal conflict 
with the improvement programme itself, since the criterion of small 0(a 2 ) corrections 
has only been touched when selecting a definite set of kinematical variables to formulate 
the respective improvement conditions within the Schrodinger functional. 

To summarize, the scaling tests in hand illustrate that also in the O(a) improved 
theory the remaining 0(a 2 ) discretization errors have to be assessed — and consequently 
can be extrapolated away reliably — by varying the lattice spacing. 

An extension of the present study to physically large volumes, where hadronic 



masses and matrix elements can be computed, is in progress [19]. 
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